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The goal of this work is to show that a ferromagnetic-like domain growth process takes place 
within the backbone of the three-dimensional ±J Edwards- Anderson (EA) spin glass model. To 
sustain this affirmation we study the heterogeneities displayed in the out-of-equilibrium dynamics 
of the model. We show that both correlation function and mean flipping time distribution present 
features that have a direct relation with spatial heterogeneities, and that they can be characterized 
by the backbone structure. In order to gain intuition we analyze the pure ferromagnetic Ising model, 
where we show the presence of dynamical heterogeneities in the mean flipping time distribution that 
are directly associated to ferromagnetic growing domains. We extend a method devised to detect 
domain walls in the Ising model to carry out a similar analysis in the three-dimensional EA spin 
glass model. This allows us to show that there exists a domain growth process within the backbone 
of this model. 



I. INTRODUCTION 



From a first approximation, glasses are characterized 
by two main features: a spatially disordered structure 
resembling the high temperature state and a slow dy- 
namics with the typical equilibration time scale exceed- 
ing the accessible laboratory time scale. It has been a 
matter of intensive study during the past decades what 
are the microscopic mechanisms behind the nature of the 
low temperature glass phase. In order to do that differ- 
ent methods and techniques have been used from both 
experimental and theoretical approaches. In particular, 
during the last years the concept of heterogeneity has 
been increasingly used to study the slow relaxation of 
glassy systems, both for structural and spin glasses. For 
structural glasses, this has led to the identification of slow 
and fast groups of particles which are spatially correlated 
and evolving with time* 1 -'^' 4 - For spin glasses, dynamical 
and spatial heterogeneities have been intensively studied, 
both from a coarse grained perspective^&£i£ where the 
dynamics is analyzed within a time or space window, or 
from single spin analysis where each spin is analyzed 
independently for a given disorder realization. 

In particular, the ground state (GS) topology of the 
three-dimensional (3D) ± J Edwards- Anderson (EA) spin 
glass modelii presents intrinsic spatial heterogeneities. 
For a given disorder realization, the multiple-degenerated 
GS configurations can be used to identify a constrained 
structure, usually referred as the "backbone" . 12 i 13 We 
have recently incorporated this spatial heterogeneous 
character of the GS into the analysis of dynamical hetero- 
geneities present in the out of equilibrium dynamic o 14 ! 15 
and also in the study of equilibrium properties of the sys- 
tem^ In Ref. we have analyzed the two-dimensional 
EA model and we have correlated the presence of slow 
and fast sets of spins with the backbone structure. Be- 
sides, for the 3D EA model we have analyzed the out-of- 
equilibrium violation of the fluctuation-dissipation the- 



orem. The results suggest that a domain growth pro- 
cess might be present within the backbone, while those 
spins outside the backbone tend to equilibrium behav- 
ior »i£ Also, by using the damage-spreading technique, we 
have results showing that ferromagnetic-like order grows 
in the 3D EA model below the pure ferromagnetic critical 
temperature T c , pointing to the presence of a Griffiths- 
like phase in the range T g < T < T c , 16 where T g is the 
glass transition temperature. Again, these results are 
suggestive of a domain growth process inside the back- 
bone structure. 

In the present work we further pursue the origin of the 
domain growth signatures found in the ±J EA model 
by using exhaustive numerical simulations. With this 
in mind we analyze the out of equilibrium dynamics in 
the 3D case while using information of the GS topol- 
ogy as a main input. By directly comparing with the 
dynamics of the pure ferromagnetic Ising model we are 
able to show strong evidence supporting the presence of 
a ferromagnetic-like domain growth process within the 
largest cluster of the backbone for T < T g . 



Outline of the paper 

Here, we outline the paper in order to guide the reader 
through the presentation of the main results. Section ITT1 
presents the model and the main characteristics of its 
GS topology. This includes the important definition of 
the backbone, allowing to classify the spins based on its 
solidary and non-solidary character. This is a useful no- 
tion to characterize spatial heterogeneities and is used 
throughout the paper. In Sec. Mil we show how these 
spatial heterogeneities are present in the two-times cor- 
relation function, establishing a direct link with typical 
finite-temperature dynamical heterogeneities. Then, the 
strong heterogeneous character of the mean flipping time 
distribution is analyzed in Sec. II VI 

Since the dynamics within the set of solidary spins 
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prompt to an analysis in terms of domain growth, we 
compare it with the well know pure ferromagnetic Ising 
model. Then, in Sec. [V] we present results for the out-of- 
equilibrium mean flipping time distribution of the Ising 
model and its correlation with the domain growth process 
which takes place. 

In Sec. IVII we describe, through a gauge transfor- 
mation, how ferromagnetic-like order can be sustained 
within the backbone structure of the 3D ± J EA model. 
It is shown that the percolative set of solidary spins can 
be thought of as a ferromagnetic-like system. All the in- 
formation is gathered in Sec. lVIIl where it is shown, after 
a comparison of the mean flipping time distribution of the 
EA model with the Ising model, that a ferromagnetic-like 
structure is growing inside the backbone of the 3D EA 
model. Section lVIIIl is devoted to the discussion and con- 
clusions. Finally, in the Appendix we characterize finite 
size effects by analyzing the out-of-equilibrium correla- 
tion function. 



II. MODEL 

We consider the 3D ± J EA spin glass model, 11 which 
consists of a set of N Ising spins crj = ±1 placed in a cu- 
bic lattice of linear dimension L, with periodic boundary 
conditions in all directions. The Hamiltonian is 



H 



(1) 



where indicates a sum over the six nearest neigh- 

bors. The coupling constants or bonds, Jij's, are inde- 
pendent random variables drawn from a bimodal distri- 
bution, i.e. Jij = ±1 with equal probability. 

We shall focus on the out-of-equilibrium dynamics of 
the model. The initial condition, t = 0, corresponds 
to a completely disordered spin configuration, that is to 
Ui = ±1 randomly chosen, which mimics a quench from 
T = oo. The working temperature is set to T < T g = 
1.12J2 Using a two-times protocol, after a waiting time 
t w we measure different time correlation functions, which 
indeed depend on the two times, t and t w . We use stan- 
dard Glauber dynamics with sequential random updates 
using a continuous time Monte Carlo algorithm. 18 

The EA spin glass model has a highly degenerate GS. 
From all the GS configurations corresponding to a single 
disorder realization it is possible to identify a constrained 
structure, which is called the rigid lattice i 12 i 13 This struc- 
ture is composed of those bonds which do not change its 
state -satisfied or frustrated- in all GS configurations. 
Using this information one can divide the spins of the 
system in two main sets: solidary spins, defined as those 
connected through the rigid lattice and thus maintaining 
their relative orientation in all GS configurations, and 
non-solidary spins, which are simply the complementary 
set. Both, the rigid lattice and the set of solidary spins 
compose the backbone of the system. 
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FIG. 1: (Color online) Correlation functions for L = 8 system 
at T = 0.6. The figure shows the full correlation function 
(black squares), the solidary (red circles) and non-solidary 
(blue triangles) contribution. The average were carried out 
over M = 10 3 samples with m = 10 thermal histories for each 
sample. Curves correspond to t w = 10 . 



The topological properties of the backbone have been 
intensively studied recentl y 19 i 20 Using a parallel temper- 
ing Monte Carlo algorithm it is possible to systematically 
arrive at GS configurations, 21 necessary to determine the 
backbone structure. We stress that the determination 
of backbone structures for each disorder realization re- 
quires to visit O(N) times the GS. In two dimensions the 
backbone does not percolate, which is consistent with the 
absence of a finite glass transition temperature. On the 
other hand, in three dimensions this constrained struc- 
ture percolates and comprises 76% of the of the spins. 
Through extensive numerical simulations we were able 
to obtain the backbone structure for M = 1000 differ- 
ent sample, i.e. different realizations of bond disorder, 
with N = 8 3 spins. This allows us to measure quantities 
that are averaged over many disorder configurations and 
thus to present results whose interpretation has a general 
physical meaning and is not restricted to a particular re- 
alization. In order to determine finite size effects we also 
compare the results obtained with L = 8 with larger sys- 
tem sizes when possible, and restrict our analysis to the 
parameter region where these effects are negligible. 



III. TWO-TIMES CORRELATION FUNCTION 

We shall first focus on the spin two-times correlation 
function defined for t > t w as 



C(t,t w ) = 



i 1 jv \ 

-5>«(£)<(iw) 
\ i=l / 



(2) 



3 



1 

0.9 
0.8 
0.7 
0.6 

0.5 

.0-4 
0.3 

0.2 



0.1 



T= 0.6 

— n— t m - 10 2 

— o — / ,= 10 3 

— A— t w = 10 4 

-v— / = 10 5 
— o— t w = 10' 

— * — t = io 7 



10 10 IO" 10 10 10 10 10 10 

/- / 




n n(b) 



10 10 10" 10 10 10 10 10 10 

/- / 



FIG. 2: (Color online) Correlation functions restricted over 
the sets of (a) solidary and (b) non-solidary spins. Curves 
corresponds to i = 8 at T = 0.6 and for different values of 
t w as indicated. The average were carried out over M — 10 3 
samples with m = 10 thermal histories for each sample. 



where af are Ising spins of sample a, (...) is an aver- 
age over m thermal histories, that is, over different ini- 
tial conditions and realizations of the thermal noise, and 
[. . . ] av indicates average over M samples. For clarity, we 
keep the a index for single disorder realizations. In the 
Appendix we discuss finite size effects and we present ev- 
idence allowing as to restrict the analysis of two-times 
quantities to systems with L = 8. 

In order to illustrate how GS information is reflected in 
the finite temperature dynamics we show here the non- 
trivial separation of the two-times correlation function. 
While the system is evolving as a whole, one can compute 
the correlation function restricting the sum over a frac- 
tion of the whole sample. Indeed, the correlation function 
restricted over the sets of solidary and non-solidary spins 



are, respectively, 



C s {t, t w ) — 



and 



C n (t, t w ) — 



(3) 



(4) 



In these definitions N" and N" are, respectively, the 
number of solidary and non-solidary spins of sample a, 
and the sums are restricted to the corresponding regions 
of each sample. Figure [T] presents the correlation func- 
tion for T = 0.6 and i w = fO 5 , comparing the full 
correlation function (black squares) with the correlation 
function restricted to the solidary (red circles) and non- 
solidary (blue triangles) spins. As discussed in Ref. [TJ 
for the two-dimensional case, the separation is non triv- 
ial. For short times the solidary spins are highly corre- 
lated presenting a very slow decay. At longer times a 
second regime with a faster decay appears. Thus, those 
spins which are highly correlated in the GS configura- 
tions, tend to maintain their correlation in time. On the 
other hand, non-solidary spins present a rapid decay of 
the correlation at short times, followed by a second faster 
decay. From this separation one can relate solidary and 
non-solidary spins to groups of spins with slow and fast 
dynamics, which is characteristic of dynamical hetero- 
geneities. Therefore, the results show that the spatial 
heterogeneity appearing in the GS configuration is cor- 
related with the dynamical heterogeneity of the out-of- 
equilibrium finite low-temperature dynamics. 

Figures^a) and (b) show, respectively, the correlation 
functions of solidary and non-solidary spins for a sys- 
tem with L = 8 at T = 0.6 and different waiting times. 
One can observe that the correlation function of solidary 
spins behaves qualitatively similar to the full correlation 
function shown in Fig. []Ja). Note also, that for solidary 
spins, at longer waiting times the initial decay is smaller. 
This means that for longer waiting times solidary spins 
are more correlated. On the contrary, the correlation 
function corresponding to non-solidary spins, Fig. E^b), 
behaves differently: the larger the values of the waiting 
time, the more pronounced is the first decay, indicating 
that non-solidary spins become less and less correlated 
with time. 

From the analysis of the correlation function given 
above, we expect that the time scale separation al- 
ready observed in the preasymptotic regime of the two- 
dimensional EA spin glass model 1 ^ will also be present 
in three dimensions. Also, it is worth stressing that in 
contrast to the two dimensional case, in 3D the back- 
bone structure percolates through the sample^ and is 
therefore potentially capable of sustaining some kind of 
phase. In the next section we analyze the time scale sep- 
aration through the analysis of the mean flipping time 
distribution. 
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FIG. 3: (Color online) Mean flipping time distribution P as 
function of lnTf, for samples of size L — 8 for At — t w — 10 7 . 
Panel (a) presents the temperature dependence of the full P. 
Panel (b) shows for T — 0.6 the two contributions, P s and 
P n , corresponding to solidary and non-solidary spins. The 
average were carried out over M = 10 3 samples, with m = 10 
thermal histories for each sample. 



IV. FLIPPING TIME DISTRIBUTIONS 

In this section we shall exploit the analysis of the 
mean flipping time distribution as an effective way to 
unveil time scale separations. In each sample we mea- 
sure the number of flips (Nf) done by every spin within 
a time window extending from f w to t, At = t — t w . 
The mean flipping time rp for a given f w and t is de- 
fined as the time window size divided by the number of 
flips: Tp(t, t w ) = At/Np-— We measure the mean flip- 
ping time distribution for each sample and then we aver- 
age the distribution over disorder realizations to obtain 
P(r F ) = [{P a (T F ))]n, v M Due to the broadness of the dis- 
tribution we use a logarithmic scale for the argument, 
i.e. P(lnTp). Figure [3Ja) shows P as function of lnrp, 



for samples of size L = 8 and different temperatures. 
We show here the particular case of At — t w — 10 7 and 
present the time evolution below. 

The distribution corresponding to T = 1.2 shows a 
single peak at small Tp, corresponding to a more homo- 
geneous high temperature situation. As temperature is 
decreased the mean flipping time distribution becomes 
broader, and eventually develops a bimodal shape, with 
a clear separation of time scales. This strong dynami- 
cal heterogeneity was shown for the first time by Ricci- 
Tersenghi and Zecchina. 2? It is also worth stressing that 
besides the clear bimodal shape, the mean flipping time 
distribution presents an internal structure characterized 
by several small and sharp peaks. These peaks, that are 
also present in Ref. |22|, follow the general temperature 
dependence of P. However, their origin is not clear, and 
will be left as an open question for future work. 

Now, using the information of the GS for L = 8, we 
show in Fig. E^b) the mean flipping time probability dis- 
tribution constrained to the sets of solidary and non- 
solidary spins, P s and P n , respectively. As was already 
observed for the two-dimensional case,— the main con- 
tribution of P s and P n are to the slow and fast peaks of 
the full P. 

Although we conclude from the study of the two-times 
correlation function that solidary spins are more corre- 
lated in time, this does not imply that they are always 
slow. Indeed, a solidary spin of a given sample does not 
always necessarily contribute to the slow right peak. This 
is also true for non-solidary spins, which can contribute to 
both, fast and slow peaks of the mean flipping time distri- 
bution functions. This is more evident if one follows the 
evolution of the mean flipping time distribution, which 
we present in Fig.|4]for T = 0.5. For clarity, we only show 
data for At = t w , the results being qualitatively similar 
for other values. Figure [4j a) shows the full P, where it 
can be observed that the bimodal shape is conserved and 
the valley between the two peaks is deeper at shorter 
times. In Figs. BJb) and (c) we show the non-solidary 
and solidary contributions, respectively, to the mean flip- 
ping time distribution. Non-solidary spins clearly present 
a second peak related to slow flips, which indicates, as 
expected, that the set of non-solidary spins is not com- 
pletely equivalent to a high temperature paramagnetic 
phase, which would be composed only of fast spins. Be- 
sides, P n slowly varies with time, the first peak and the 
valley between peaks being higher for longer times. In 
Fig. [4|c) one can observe that P s also presents two peaks 
at short times. However, whilst the second peak does not 
change with time, the first peak presents a drastic fall. 
This is a strong and important difference between the 
relaxational behavior of solidary and non-solidary spins. 
A clue to understand this behavior comes again from the 
information contained in the GS configurations. Before 
going into more detail we will analyze a canonical model 
in statistical physics: the Ising model. This will allow us 
to gain intuition and also to determine the key observ- 
ables necessary for a quantitative analysis. 
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FIG. 4: (Color online) The mean flipping time probability distribution at T — 0.5 for different times and waiting times as 
indicated, keeping At = t w . Different panels correspond to: (a) all spins, P, (b) non-solidary spins, P n , and (c) solidary spins, 
P B . The average were carried out over M = 10 3 and m = 10 for L = 8. 
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FIG. 5: (Color online) The mean flipping time distribution 
of the 3D Ising model at T = 2.0 for different values of t w 
with At = t w . The average were carried out over m = 5 x 10 4 
thermal histories for L = 10. 



V. THE ISING MODEL 

In this section we use the same protocol defined in the 
previous section to analyze the dynamics of the pure 3D 
ferromagnetic Ising model, where all bonds are identical, 



-J. In particular, we simulate a quench from 



T = oo to a working temperature T < T c = 4. 515. 23 

In Fig. [5] we present the time evolution of P(lnrp) 
when the system is quenched to T — 2.0. We use here 
L = 10 in order to have a system size comparable to the 
one used for the EA model. We also used L = 50 and ob- 
tained the same qualitative behaviour. For short times 
the distribution presents a clear bimodal shape. How- 
ever, for increasing i w , the behavior is different from the 
one observed in the full mean flipping time distribution 
of the EA model, Fig. @fa). In particular, the peak on 



the left rapidly falls with increasing waiting times. This 
behavior strongly resembles the one observed for solidary 
spins in Fig. Hfc). 

The evolution of the Ising model below its critical tem- 
perature serves as a canonical example of growing ferro- 
magnetic order. In order to advance further in under- 
standing the behavior of the mean flipping time distri- 
bution we will present a quantitative characterization of 
the domain wall dynamics. This can be done by imple- 
menting the algorithm of Hinrichsen and Antoni, 24 which 
permits to identify for each of the thermal histories which 
spins flip due to a domain wall. The main idea of the al- 
gorithm is to simulate simultaneously three replicas of 
a system with different initial conditions (disordered, all 
spins up, and all spins down) with exactly the same ther- 
mal noise, i.e. using the same random number sequence 
for the three replicas. The spins in the replicas that be- 
gin with an ordered configuration will only flip due to 
thermal noise. As a consequence, a simultaneous flip in 
all three replicas is considered as due to thermal noise. 
On the other hand, if a spin does not flip simultaneously 
in all three replicas it is due to a domain wall. 

We can thus determine Nw, the number of time steps 
that a spin spends in a domain wall within the time win- 
dow At. Thus we can compute the distribution function 
of the mean time rw = At/iVw that a spin spends in a 
domain wall, i.e. the mean domain- wall-time distribution 
QQhtw). Note that this is also an out-of-equilibrium 
quantity depending on t and t w . Figure [5] shows the evo- 
lution of Q(lnrw) at different times, and also a compari- 
son with P(lnrp). As can be observed, Q is also a highly 
heterogeneous function. For short time scales the left 
peak of both distributions coincide. For increasing times 
the valley in Q becomes broader, and eventually the dis- 
tribution becomes almost homogeneous beyond the first 
peak, which is always present. This behavior is somehow 
mirrored by P, that, as already described above, presents 
a decay in the fast peak. 

The physical interpretation of the behaviors of P and 
Q is intimately related to the domain growth process in 
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FIG. 6: (Color online) Mean flipping time distribution P vs. 
hiTF (open square) and mean domain- wall-time distribution 
Q vs. lnrw (solid circle), of the 3D Ising model at T = 2.0 for 
different values of t w with At = i w . The average were carried 
out over M = 5 x 10 4 for L = 10. 



the Ising model. For low temperatures most flips will be 
due to domain walls, while flips inside a domain will take 
place with a very low probability. For short times the 
domains are very small, and the dynamics is fast in the 
sense that there are many flips due to domain walls dy- 
namics, and both distributions coincide. For increasing 
time, domains grow and the fraction of spins in domain 
walls decay. This is reflected in the behavior of both P 
and Q for small values of, respectively, rp and rw, and 
establishes a relation between dynamical heterogeneities 
and domain growth in the Ising model. Therefore, the 
strong correlation between P and Q observed in Fig. [5] 
indicates that the characteristic decay with time of the 
first peak is a hallmark of a domain growth process. 

The first question that comes immediately to mind is 
if it is possible to relate some kind of phase ordering 
to the dynamical heterogeneities observed in the 3D EA 
model as we just did for the Ising model. In the following 
sections we will tackle this issue. 



VI. GAUGE TRANSFORMATION 

In this section we shall present evidence that supports 
the idea that the backbone can sustain a ferromagnetic- 
like ordered structure.— Recall that in Sec.|H]we pointed 
out that in three dimensions the backbone percolates. 
Then, here we will focus our interest by restricting the 
following analysis to the largest cluster of solidary spins. 
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FIG. 7: (Color online) Two-dimensional example of the gauge 
transformation, (a) The GS configuration of solidary spins: 
black (white) dots represents up (down) spins, single (dou- 
ble) lines correspond to ferromagnetic (antiferromagnetic) 
bonds, and red (blue) color corresponds to frustrated (sat- 
isfied) bonds, (b) After the gauge transformation all spins 
are up, all ferromagnetic bonds are not frustrated and all an- 
tiferromagnetic bonds. 



Thus, we consider the Hamiltonian 

(*,i)e{«} 



(5) 



where the sum has been restricted to this percolating set 
of solidary spins. Then, we consider the new spin vari- 
ables denned through Sj = ofoi, where {of} represent 
one of the two possible values of the spins variables in 
all the GS configurations. Now, by replacing the bonds 
by = Jij<j®<Tj, and since (of) 2 = 1 one arrives at the 
Hamiltonian 



(*.j)£{s} 



J > 



(6) 



which is indeed equivalent to H s . Therefore this gauge 
transformation leaves the Hamiltonian invariant. It is 
clear that in the GS the new spin variables are all = 
+1, or equivalently the configuration with the opposite 
sign. 

As an illustrative example, consider the set of solidary 
spins in Fig. [TJa) , which corresponds to one of the two 
possible GS configurations of the largest cluster of the 
backbone in a particular realization of a two-dimensional 
system. In the figure, single (double) lines correspond 
to ferromagnetic (antiferromagnetic) bonds from the bi- 
modal distribution; red or blue color indicates if the in- 
teraction is frustrated or not, respectively. Figure [7£b) 
represents the new set of spins {s°} with the correspond- 
ing bonds J-j. Notice that after the transformation all 
spins have the same sign. Another property of the 
gauge transformation is that it preserves the frustration 
in the system. Note that all the frustrated (red) bonds 
in Fig. [7] are always frustrated. Indeed, in three dimen- 
sions approximately only a 10% fraction of the bonds 
linking spins within the backbone are frustrated. In this 
way we have obtained, after the gauge transformation, 
a ferromagnetic-like system with a 10% fraction of frus- 
trated antiferromagnetic bonds. Since 0.1 < x c — 0.222, 
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FIG. 8: (Color online) The distributions P s vs. hiTF (open 
square) and Q a vs. lnrw (solid circle) corresponding only to 
the subset of solidary spins of the EA model. Data correspond 
to T — 0.5 for different values of i w with At = t w . The average 
were carried out over M = 10 3 and m = 10 for L = 8. 



the limit of concentration of antiferromagnetic bonds 
in the ferromagnetic phase, 25 and since the backbone 
structure percolates in three dimensions,^™ we expect 
that the set of solidary spins is a good candidate for 
ferromagnetic-like order. In particular, we expect that 
the qualitative behavior of the evolution of P and Q for 
solidary spins will somehow resemble the ones observed 
in the Ising model. 



and dictates how the systems can be decomposed in sol- 
idary and non-solidary spins. Since, by definition, sol- 
idary spins maintain their relative orientation in all the 
ground state configurations we can choose just one GS 
configuration -and the one with all the spins reversed- 
to compare with the replica that has a random initial 
condition. In this way the algorithm will permit us to 
identify domain walls within the backbone, which is the 
structure with a simple symmetry break. On the other 
hand, we will not be able to obtain any relevant informa- 
tion outside the backbone. Once the domain walls have 
been identified one can compute the mean domain-wall- 
time distribution, Q~2- 

In Fig. [8] we present the mean flipping time distribu- 
tion for solidary spins P s vs. mrp (open square), to- 
gether with the mean domain-wall-time distribution for 
solidary spins Q s vs. lnrw (solid circle) for T = 0.5. 
As can be observed, both distributions are highly corre- 
lated, presenting similar qualitative behaviors. The peak 
on the right does not seem to depend on the waiting 
time. On the other hand, for short waiting times the 
distributions present a peak on the left, which decreases 
with increasing waiting time, and eventually disappears. 
Therefore, the mean flipping time and the mean domain- 
wall-time distributions for solidary spins are highly cor- 
related and presents a first peak which displays a fast 
decay with increasing waiting time. As was observed for 
the Ising model in Sec. [V] this behavior is intimately re- 
lated to a domain growth proccess. Thus, the comparison 
with the Ising model and also the small bond frustration 
present on the backbone suggest a plausible scenario for 
ferromagnetic-like growing order in the 3D ± J EA model. 



VIII. DISCUSSIONS AND CONCLUSION 



VII. DOMAIN GROWTH WITHIN THE 
BACKBONE 

In this section we shall present further numerical evi- 
dence for a growing ferromagnetic-like process in the EA 
model. In Sec. |V] we analyzed the mean flipping time 
distribution for the Ising model. We observed that this 
distribution presents a bimodal shape, and we were able 
to establish a direct relation between the time that spins 
spend in a domain wall and the changes in shape that 
takes place in the mean flipping time distribution as a 
domain growth process takes place. 

In order to establish a direct comparison with the re- 
sults of the previous sections we will present a gener- 
alization of the algorithm for domain wall detection of 
Hinrichscn and Antonij^l originally devised for the Ising 
model. Since the ±J EA model has a highly degener- 
ate GS, the perspective of simulating the system and 
replicas in all the GS configurations seems a daunting 
task. However, as has been stressed through the article, 
the ground state of the 3D EA model is very peculiar, 



Throughout the present work, based on the separation 
in solidary and non-solidary spins given by the GS topol- 
ogy, we have presented evidence supporting a picture of 
growing ferromagnetic-like order within the backbone of 
the 3D ± J EA model. 

By studying the two-times correlation function we 
showed that the solidary/non-solidary spin separation is 
not trivial in the sense that it is able to give meaning- 
ful dynamical information. While solidary spins maintain 
their correlation in time, non-solidary spins present a fast 
initial decay of the correlation function. 

On the other hand, the mean flipping time distribu- 
tion shows that solidary (non-solidary) spins mainly con- 
tribute to the slow (fast) peak. Besides, although it can 
be shown that a given spin contributes to both peaks, the 
time evolution of the mean flipping time distribution for 
solidary and non-solidary spins behaves qualitatively dif- 
ferent. As a main characteristic, the left peak of the sol- 
idary spins P s rapidly decreases with increasing waiting 
time. In order to correlate this information with ferro- 
magnetic order we have also analyzed the pure ferromag- 
netic Ising model using the same protocol. Thus, we were 



8 



able to associate the heterogeneous behavior of P(Iiitf) 
with the mean domain- wall-time distribution <5(lnT\v)- 
The results obtained are the expected for a coarsening 
system like the Ising model and give us a quantitative 
tool to analyze the mean flipping time distribution in 
the EA model. 

Seeking for ferromagnetic order in the 3D ± J EA spin 
glass model, we showed that using a gauge transforma- 
tion the solidary spins can be transformed into a ferro- 
magnetic system with a small fraction of frustrated anti- 
ferromagnetic bonds. With this in mind we presented a 
generalization of the algorithm for domain wall detection 
of Hinrichsen and Antoni. This allowed us to present a 
quantitative characterization of P s (tf) and Q s (tw) and, 
based in the comparison with the Ising model, we con- 
clude that ferromagnetic-like order is growing within the 
backbone of the model. 

All these results, together with other recent works that 
show the relevance of the heterogeneous character of the 
GS topology ; 14 ' 15 i 16 ' 26 point to an alternative picture for 
the nature of the spin glass phase. Within this pic- 
ture, ferromagnetic-like order grows within a constrained 
structure of the system. It was precisely with this idea in 
mind that we were able to present a generalization of the 
method proposed by Hinrichsen and Antoni, and thus 
characterize a ferromagnetic-like domain growth process. 

Finally, we stress the non-trivial character of the re- 
sults by comparing with the Gaussian spin glass model 
that has only one simple degenerated GS. Note that even 
when a gauge transformation can lead to a disorder fer- 
romagnetic system, the generalization of the domain wall 
detection method is not direct due to the high frustration 
of the whole system. In order to actually obtain relevant 
physical information it would be thus necessary to deter- 
mine some kind of low frustrated constrained structure. 
Concluding, the results presented in this work highlight 
the importance of the determination of the backbone 
structure in spin glass models. 
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APPENDIX: FINITE SIZE EFFECTS 

We shall show within this appendix that finite size 
effects are not strong in the time window we are in- 
terested in through the paper. For this purpose, al- 
though one can find the temperature and time depen- 
dence of the two-times correlation function in previous 
works, 27,28 in the following we use the out-of-equilibrium 
correlation function to carefully compare different system 




FIG. 9: (Color online) Two-times correlation function for 
L = 8 (open symbols) and L = 20 (solid symbols) systems 
at (a) T — 0.6 for different values of t w , and at (b) £ w = 10 5 
for different temperatures as indicated. The averages were 
carried out over M = 10 3 (L = 8) and M = 10 2 (L = 20) 
samples, with m = 10 thermal histories for each sample. 



sizes. Figure [Ha) shows C(t,t w ) for a fixed temperature 
T = 0.6 < T g and different waiting times t w . Data are 
shown for two system sizes, L = 8 and L = 20. Up 
to delay-times of the order of t — t w ~ 10 5 the curves 
overlap. For larger delay-times small deviations are ob- 
served due to finite size effects. However, we stress that 
the same qualitative behavior is observed for both sys- 
tem sizes. The correlation function for the same system 
sizes is also considered in Fig. |5fb), but as a function of 
temperature for fixed t w = 10 5 . It is observed that finite 
size effects are stronger at higher temperatures. However, 
these are not significant for temperatures below T = 0.6 
and times such that At = i w = 10 5 . For small tempera- 
tures and larger times there are small deviations between 
the L = 20 and L = 8 curves, which nevertheless show a 
good qualitative agreement. 
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FIG. 10: (Color online) Correlation function for L = 8 sys- 
tems at T = 0.6, calculated with Eq. ((2)1 (open symbols) and 
Eq. (|A.1[) (solid symbols). The average were carried out over 
M = 10 3 . 

A different approach we use to quantify finite size ef- 
fect comes from the information obtained with the sepa- 
ration of the correlation function in its solidary and non- 
solidary contributions, as explained in Sec. Defining 
the mean fractions of solidary and non-solidary spins as, 
respectively, / s = [N?/N] &v and / n = [N^/N] av with 



fs + fa = 1: then in the thermodynamic limit it should 
hold that 

C(t,0 - f s C s (t,t w ) + f n C n (t,t w ). (A.l) 

This would be strictly valid in the thermodynamic limit. 
If finite size effects are important, the mean fractions 
would present important fluctuations and Eq. ([2]) could 
not be split in the solidary and non-solidary contribu- 
tions. We show in Fig. \TU\ that, although the mean 
fractions do fluctuate^ the total correlation function is 
very well approximated by Eq. (|A.1[) . even for L = 8. 
Data correspond to the full correlation function com- 
puted using Eq. @ (open symbols) and the one com- 
puted from Eq. (|A.1[) (solid symbols). Notice that the 
different data sets are almost indistinguishable, even for 
the small system size we are using. In this case, averages 
were taken over M — 10 3 independent samples which 
yields / s = 0.75883 and /„ = 0.24117. 

Summarizing, in this appendix we have performed a fi- 
nite size study by comparing two system sizes, L = 8 and 
L = 20. We have shown that for the low temperatures 
and for the time scales we are interested in, physically 
relevant information can be extracted from a system of 
linear size L = 8. Thus, to analyze two-times dependent 
observables of the 3D ± J EA model we used this system 
size, where we have determined the backbone structure 
for an extensive number (M = 1000) of disorder realiza- 
tions. 
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We have compared the results obtained using one GS changes, 
configuration with the one obtained by averaging over 
many GS configurations and we checked that there are not 



